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STUDY OF THE STABILITY AND DRAG AT MACH NUMBERS 


FROM 4-5 TO 13-5 OF A CONICAL 

VENUS -ENTRY BODY 

By Peter F. Intrieri 
Ames Research Center 


SUMMARY 


An experimental study has been conducted to determine the static and 
dynamic stability and drag characteristics of a conical -nosed planetary-entry 
vehicle in a gas mixture composed of 9”P ercen ^ carbon dioxide and 91 ~percent 
nitrogen and in air at Mach numbers of 4.5* 9*0, and 13*5 and Reynolds numbers 
from 0.4 to 1.2 million. It was determined that for the assumed center-of- 
gravity location, the configuration is statically and dynamically stable at 
all Mach numbers and throughout the angle -of -attack range of the investigation 
and that it exhibits the same static stability, dynamic stability, and drag in 
air as it does in the nitrogen --carbon -dioxide mixture . The static stability 
is invariant with Mach number and nonlinear with angle of attack. The non- 
linear variation of pitching moment with angle of attack was closely approx- 
imated by a two -term power series of the resultant angle of attack. The 
dynamic stability is invariant with Mach number and angle of attack. The drag 
coefficient of the configuration is essentially invariant with Mach number and 
increases with increasing angle of attack. The static stability and drag were 
closely predicted by Newtonian theory. 

Strong experimental evidence is presented which shows that the configura- 
tion is statically stable only about a nose -forward trim attitude of zero 
angle of attack and dynamically stable about this unique stable attitude for 
all amplitudes of oscillation. 

Calculations which use the measured stability results to describe the 
oscillatory behavior of a full-scale vehicle flying an example entry trajec- 
tory through a model of the Venus atmosphere indicate that the stability char- 
acteristics of the present configuration should be adequate for orienting the 
vehicle properly during a Venus entry. For all initial orientations of the 
vehicle (nose forward or base forward) the pitching motions should converge to 
a very small fraction of the amplitude at entry. 


INTRODUCTION 


In the design of vehicles intended to enter planetary atmospheres, two 
important goals are to provide for minimum aerodynamic heating and for ade- 
quate aerodynamic stability of the vehicle during entry. For entries into the 
Earth 1 s atmosphere at velocities well in excess of escape velocity. 
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reference 1 has shown that the total heating will be less for sharp entry 
bodies, such as pointed cones, than for blunt bodies- Radiative heating which 
becomes the dominant mode at these speeds depends on the velocity normal to 
the bow -shock wave rather than on the free -stream velocity; therefore, the 
combined radiative and convective heating will be drastically reduced if the 
entry bodies have highly swept bow-shock waves. Reference 2 indicates that in 
the nitrogen --carbon -dioxide atmospheres of Venus and Mars, the radiative 
intensities will be significantly greater than in air. The advantages of 
using conical entry bodies to reduce total heating in these atmospheres are 
shown by the analysis of reference 3 • 

Having established the desirability of conical shapes from heating 
considerations, it must be ascertained, for a passive system, that the stabil- 
ity of such vehicles flying in the atmospheres will be adequate to 

insure proper orientation of its heat shield during the entry, and, in partic- 
ular, that the vehicle* s dynamic stability will be sufficient to prevent 
divergent oscillations which could result in a mission failure. Also, since 
the constituents of the planetary atmospheres through which these vehicles 
must descend are known only within certain limits (see, e.g., ref. 4), it is 
important to determine whether the aerodynamic characteristics of the vehicle 
are sensitive to changes in gas composition and whether these characteristics 
can be predicted satisfactorily. 

The static and dynamic stability and drag characteristics of a conical 
configuration were investigated in the prototype of the Ames Hypervelocity 
Free -Flight Facility at nominal Mach numbers of 4.5, and 13*5 and at 

Reynolds numbers of 0.4, 0.8, and 1.2 million, respectively, based on free- 
stream conditions and model diameter. The models were tested in a gas mixture 
composed of 9~percent carbon dioxide and 91 - percent nitrogen (by volume) and 
also in air. Two exploratory flights were also made with the models launched 
backward to determine whether they would remain flying backward or, as 
expected from calculations, would begin righting themselves to a nose -forward 
attitude. Theoretical estimates of static stability and drag were made and 
compared with the experimental results. 


SYMBOLS 


A reference area, maximum body cross-sectional area, meters 2 

C-l arbitrary constant in equation (10) 

Cp drag coefficient, — 9 dimensionless 

( 4xA 

Cp Q drag coefficient at zero angle of attack 

Ct lift coefficient, — , dimensionless 

Ct lift-curve slope , per radian 
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pitching -moment coefficient, P^ c ^ in g moment dimensionless 

pitching -moment -curve slope , per radian 

5Cm SC™ 

damping -in -pitch derivative, — j— + —— t — j— , per radian 

o( qd/V ) o(ad/V) 

reference diameter, maximum body diameter, cm 

average gravitational acceleration within the atmosphere of a 
planet, m/sec 2 

altitude above planet surface, km 

moment of inertia about a transverse axis through the center of 
gravity, kg-m 2 

constant in equation ( l) 

dynamic -stability parameter for variable density (eq. (ll)), 
dimens ionle s s 

constants in equation (2), deg 

mass of model, gm 

Mach number, dimensionless 

re storing -moment coefficients, defined by equation (Al) 

re storing -moment coefficients, defined by equation (A4) 

integer denoting number of complete tumbles of vehicle 

roll parameter, roll rate/velocity, radians/m 

static pressure, newtons/m 2 

angular pitching velocity, radians/sec 

free -stream dynamic pressure, newtons/m 2 

gas constant for atmospheric gas mixture, , m 2 /sec 2 °K 

Reynolds number based on free -stream gas properties and maximum 
diamete r , dimens ionle s s 


time 
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T gas temperature , °K 

u horizontal component of flight velocity, m/sec 

V velocity along flight path, m/sec 

x distance along flight path, m 

x„~ axial distance from model nose to center -of -gravity position, cm 

y horizontal coordinate normal to the flight path, m 

z coordinate normal to the flight path and y axis , m 

a angle of attack (angle between model axis and resultant wind direc- 

tion projected onto the vertical plane) , deg 

3 angle of sideslip (angle between model axis and resultant wind 

direction projected onto the horizontal plane), deg 

3^ reciprocal of atmospheric scale height, ^ , per m 

7 flight -path angle (referenced to the local horizontal) , deg 

e dependent variable, - ^ + nit, deg 

€p value of e at the first peak of oscillatory motion 

damping exponents in equation (2) , m -1 
0 ,i[r attitude coordinates of the model relative to Earth -fixed axis, deg 

A wave length of pitching oscillation, m/cycle 

| dynamic -stability parameter for constant altitude (eq. (5)), dimen- 

sionless 

p atmospheric density, kg/m 3 

a resultant angle of attack, tan -1 0 tan 2 a + tan 2 3 , deg, or transverse 

Iy ( (j 

radius of gyration, -g- , when used in parameter l in equa- 
tions (5), (8), and (ll) and in table I, m 


°m 

maximum resultant angle of 

attack, deg 


a o 

minimum resultant angle of 

attack, deg 




i 

n* 

Q 

ro 

a rms 

root -mean -square re suit ant 

angle of attack, / 

x dx, deg 

GO 

frequency , radians / sec 
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^ 2 .^ 2 . rate of rotation of complex vectors which generate the model pitching 
motion (eq. (2)), radians/m 

( ) first derivative with respect to time 

( ) second derivative with respect to time 

Subscripts 

eg center of gravity 

env envelope 

i initial condition 

max maximum 

min minimum 

o conditions at surface of planet 

oo free -stream conditions 


CONF IG-URATION 


The configuration selected had a 3 O 0 half -angle cone forebody. This 
forebody , as is shown in references 1 and 3 ; is close to the optimum cone 
angle for minimizing total heating at velocities well above escape velocity. 
The afterbody geometry was selected by the requirement that the configuration 
be stable only with the 30° half -angle cone pointed forward. This requirement 
was fulfilled by the addition of a short cylinder and a rearward pointing cone 
to the base of the forebody. The cylinder length of 0.125 model diameter and 
the rearward pointing cone half -angle of 48° were selected to provide a coin- 
cident center -of -pressure location for both base-forward and nose-forward 
flight of the vehicle and also to provide useful volume. Figure l(a) is a 
sketch of a model of this configuration, giving pertinent nominal dimensions. 


DESCRIPTION OF TESTS 


Test Technique and Test Conditions 

The tests were performed in the prototype of the Ames Hypervelocity Free- 
Flight Facility by launching models from a deformable -piston, light -gas gun of 
12.7 mm bore (ref. 5 ) into a quiescent mixture composed of 9 tpercent carbon 
dioxide and 91 “Pe rcen 'k nitrogen (by volume), and into still air at ambient 
temperature. Nominal model velocities of 1*52, 3-05, and 4.72 km/sec were 
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obtained corresponding to nominal Mach numbers of 4.5, 9.0, and 13.5, respec- 
tively. In order to obtain adequate definition of the motions of the models 
through the test section, with the given spacing of shadowgraph stations (see 
below), the ambient test-section pressure was set at 4xi0 4 newtons/m 2 (300 mm 
of mercury), to give the desired wave length of oscillation. Corresponding 
Reynolds numbers are shown in table I, which lists the average values of Mach 
number, Reynolds number, and test gas for each flight. Gas samples taken 
immediately prior to each firing in the N 2 -C0 2 mixture were analyzed for chem- 
ical content on a mass spectrograph and were found to contain the desired con- 
centration of C0 2 and N 2 within 1 percent, with only trace amounts of other 
constituents . 

The trajectory of the model through the test section was recorded by 11 
spark shadowgraph stations located at 1.22 -meter intervals. Side and top 
views of the models were recorded at each station. The shadowgraphs contained 
images of reference wires from which x,y,z ,0, and \|r coordinates were read. 
The linear coordinates were measured to within O.OO76 cm, and the angles to 
within 0.25°« The orientation angles 9 and ^ were read relative to Earth- 
fixed axes. No corrections were made for the angle between the resultant wind 
direction and Earth -fixed axes to give values of a and £ since, for these 
tests, these corrections were within the reading accuracy of the angles 9 
and Time of model flight between stations was recorded with electronic 
chronographs to within 0 .05 microsecond (0.02 percent of the time of flight 
between stations of a model flying at 4.72 km/sec). 


Models and Sabots 

Sketches of the two types of models used in the present investigation, 
showing pertinent nominal dimensions, are presented in figure 1. The after- 
sections of the models shown in figure l(a) were machined from 7075 _ T6 alu- 
minum; the front sections of the models were machined from a tungsten-iron 
alloy (Mallory 3^0°) to obtain a center -of -gravity location of 0.72 d from 
the nose. These models were used in the tests at Mach numbers of 4.5 and 9 *0 . 
However, at the highest Mach number (M = 13-5); the heating rates encountered 
were sufficiently high to cause the models to burn in flight. It was found 
that a thin plastic coating over the models provided a successful heat shield. 
The coating consisted of a polymer of vinylidene fluoride sprayed onto the 
model surface to a thickness of about 0.04 cm. The models were remachined and 
trimmed to give a final coat thickness of about 0 .019 cm. Figure l(b) is a 
sketch of the coated models showing pertinent nominal dimensions. The after- 
sections of these plastic -coated models were machined from a titanium alloy 
since the heat -curing of the plastic was found to anneal the aluminum after- 
section so that it deformed severely during the launch. The strength of the 
titanium alloy was not affected by the curing. Since titanium is heavier than 
aluminum, the size of the Mallory front section was adjusted (see fig. l(b)) 
so that the center -of -gravity location of these models, prior to being coated, 
was the same as that for the basic models - 0.72 d from the nose. The amount, 
and consequently the weight distribution, of the plastic used in the coating 
was small; therefore, the location of the center of gravity relative to any 
point on the metal model did not change. However, the distance of the center 
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of gravity to the new plastic nose of the coated model was increased by the 
length of plastic deposited at the nose - approximately 0 .O38 cm. Examination 
of the shadowgraph pictures (fig. 2 ) of two of these models in flight at Mach 
number 13-5 shows that although ablation of the plastic occurred during the 
flights, any change in shape , particularly blunting of the tip, was insignifi- 
cant . The position of the center of gravity for all the models was measured 
within ± 0.0025 cm accuracy. The dimensions of the models varied only slightly 
from those shown in figure 1 . Some of the measured physical characteristics 
of each model are listed in table I. 

Photographs of the models and sabots are presented in figure 3. Three 
different sabot types were used in launching the models: a straight sabot 

(fig. 3( a )) for flights at angles of attack below 10°, a canted sabot 
(fig. 3(b) ) for flights at angles of attack as high as ^5° , and a sabot for 
launching the models in the base -forward attitude (fig. 3(<i))- The canted 
sabot was not always successful in producing high angle -of -attack flights and, 
as can be seen in table I, resulted in a large number of low angle -of -attack 
flights in the N 2 -C 0 2 mixture . Tests using the canted sabot to launch the 
coated models at a Mach number of 13 .5 were unsuccessful because the models 
broke up inside the gun. The straight sabot was successful in obtaining low- 
angle flights at this Mach number. The sabots were machined from Lexan 
polycarbonate plastic . 


REDUCTION OF DATA 


Drag 

The determination of drag coefficient from the time -distance data was 
based on the procedure described in reference 6, in which a constant drag 
coefficient is assumed. A procedure applicable to cases where the drag coef- 
ficient varies with angle of attack is presented in reference 7 * It is shown 
in this reference that if the drag coefficient varies with the square of the 
local resultant angle of attack, according to the relation 

C D “ C D 0 + kcj2 (l) 

the drag coefficient obtained by the method of reference 6, under certain 
additional constraints, is the drag coefficient that would be obtained in 
steady flight at a resultant angle of attack equal to the root -me an -square 
resultant angle of attack, cr rms . Accordingly, the present results are cor- 
related with a rms . 


Stability Derivatives 

The pitching and yawing motions of the models during free flight were 
analyzed to determine the stability derivatives. Examples of the types of 
motions encountered in the present tests, as viewed in the a-|3 plane, are 
shown in figure 4 . The angles of attack and sideslip measured from the 
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shadowgraphs at each station are indicated by the symbols . The curves show 
the theoretical motions which best fit the experimental data and were computed 
by a method which will be discussed later in this section. Since the models 
are aerodynamic ally symmetric, the angular displacement of the model, at any 
instant, can be represented also by the resultant angle of attack, o', whose 
orthogonal components are the angles a and 0. It can be seen in figure 4 
that, in general, the data show precessing elliptical motions, and that the 
angle range through which the models oscillate differs for each flight. Also, 
it should be noted that the models in the tests at Mach numbers of ^.5 and 
9.0 underwent two complete cycles of oscillation (figs, 4(a), (b), (c), and 
(d)), whereas the models in the tests at Mach number of 13-5* because they 
were heavier (see discussion under Models and Sabots), underwent only one and 
one-half cycles of oscillation (figs. 4(e) and (f)). Although one and one-half 
cycles of motion are usually sufficient to define the desired stability deriv- 
atives, more cycles of well-defined motion will generally produce better data. 

Stability derivatives were obtained from analysis of the pitching and 
yawing motions of the models. For each flight, the measured variation of a 
and 0 with x was represented by the following equation: 


0 + ia = Kje 


(■ni +iUJ i) x 


(r] 2 -iw )x 


+ K 2 e N '' 2 + Kse ipx 


( 2 ) 


Equation ( 2 ) is the solution of the linear differential equation of motion, as 
given in reference 8, and includes effects of model spin and trim angle on the 
motion. Some of the basic assumptions used in the development of this equa- 
tion are: axially symmetric configuration, linear variations of force and 

moment with angle of attack, small angular displacements, and small angles of 
trim. Equation ( 2 ) programed for machine computation (ref. 9) was used to 
select optimum values of all the constants by an iterative process of differ- 
ential corrections. The curves shown in figure 4 were obtained by fitting 
equation (2) to the experimental data. The closeness of the computed curves 
to the experimental data is a measure of the reliability of the stability 
results. The fitted curves for all the flights analyzed in this investigation 
agreed with the measured angles within the measuring accuracy. 


The pitching -moment -curve slope, C^, ’was computed from the wave length 
of oscillation by means of 


where 




8jt2l y 

A 2 pAd 


X = 


2jt 


(3) 

( 4 ) 


The dynamic -stability parameter. 


£, defined as 


2 


(5) 


6 - C D - %* + + C m d \ a 


was determined from the constants % and r\ 2 by means of the relation 


pA t 

li + ^ * S 6 


(6) 


It has been shown in references 9 and. 10 that | , in the form shown in equa- 
tion ( 5 ) , is a convenient parameter for describing the dynamic stability of an 
unpowered vehicle in free flight at constant altitude. The values of £, pre- 
sented in this report^ were obtained from equations ( 2 ) and ( 6 ) which assume a 
linear system over the angle -of -attack range covered during any one flight. 
Each value of f- , therefore; is the dynamic -stability parameter of an equiva- 
lent linear system whose amplitude of oscillation would grow or diminish in 
the same way as that experienced by the model. 


RESULTS AND DISCUSSION 


The experimental results of this investigation define the drag; static - 
stability, and dynamic- stability characteristics of the configuration in the 
angle -of -attack range from 0 ° to 40 °; at Mach numbers of 4.5* 9 -0 ; and 13 -5 j 
in a gas mixture composed of 9 “Percent carbon dioxide and 91“Percent nitrogen 
and also in air. The measured values of Cp; and |; obtained from anal- 

ysis of 65 separate model flights; are summarized in table I. These coeffi- 
cients are based on the maximum model diameter and frontal area. Values of 
°m> °o> an d- a rms presented in table I indicate the angle -of -attack range 
through which each model oscillated during the flight. Full-scale enthalpy 
values were duplicated in these tests and; as stated earlier; at the highest 
Mach number (13-5)., the heating rates were sufficiently high to cause ablation 
of the protective plastic coating on the models. These models survived the 
flights without significant change in surface shape. The effects of gas dis- 
sociation and surface ablation cannot be isolated but are implicit in the 
experimental results. Shadowgraphs; typical of those obtained in the present 
tests; are presented in figures 5 and. 6 to illustrate gross features of the 
flow fields; particularly the bow-shock wave. Theoretical estimates of static 
stability and drag made using conical flow theory (ref. 11 ) ; and the equations 
developed in reference 12; based on Newtonian impact theory; are compared with 
the experimental results. 


Drag Characteristics 

The measured values of drag coefficient presented as a function of the 
root -mean-square resultant angle of attack in figure 7 show that the drag 
coefficient increases approximately 30 percent as the angle of attack is 
increased from 0° to 28°. These data also show that in this angle -of -attack 
range the configuration exhibits the same drag characteristics in air as it 
does in the carbon dioxide and nitrogen mix±ure . The drag coefficient is seen 
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to decrease slightly (approximately 8 percent) as the Mach number is increased 
from 4.5 to 9*0 but is seen to remain constant; within the scatter of the data; 
for a further increase in Mach number to 13 .5 • Comparison of the experimental 
results with values calculated for zero angle of attack using conical -flow 
theory shows that although the theory is in excellent agreement with the exper- 
imental data at Mach numbers 9-0 and 13*5; and predicts the correct trend with 
Mach number; the theory underestimates the drag at Mach number 4-5 by about 
6.5 percent. Estimates were made of the base-pressure contribution to the 
drag for all three Mach numbers using a Prandtl-Meyer expansion around the 
base of the model and flow deflection angles measured from shadowgraphs simi- 
lar to those presented in figures 5 and 6 . These estimates if applied would 
raise the theoretical values at Mach numbers 4.5, 9-0; and 13-5 by about 6.5; 
1 . 5 ; and 1 percent of the measured values, respectively, and would result in 
excellent agreement with the experimental results . The effect of angle of 
attack on drag is well predicted by Newtonian flow theory. 


Static -Stability Characteristics 

The experimental values of the equivalent linear pitching -moment -curve 
slope, C^, are presented in figure 8 as a function of the maximum resultant 
angle of attack, a m . For a nonlinear system the stability results, unlike the 
drag results, do not correlate with a rras so are presented as functions of 
a m , which is convenient for the further analysis of the data by nonlinear 
methods. The data have been corrected for the variations in center -of -gravity 
location shown in table I to a common center -of -gravity position, 0.72 d from 
the nose. With the exception of flight number 853 (see table I) , these cor- 
rections to the data of Mach number of 4.5 and 9*0 were small, usually within 
4 percent (the measuring accuracy of center -of -gravity position) . However, 
the corrections to the data of Mach number 13*5 were on the order of 20 per- 
cent, since, as discussed earlier, at this Mach number the models required a 
protective plastic coating which caused their centers of gravity to be at 
0.733 cL from the nose (see fig. 1 (b)). 

The data, in figure 8 , show that the configuration is statically stable 
throughout the angle -of -attack range investigated for each Mach number. These 
data also show that the static stability is insensitive to changes in Mach 
number from 4.5 to 9 *0 • Comparison of these data with the data at Mach number 
I 3.5 shows that the combined influence of further increasing the Mach number 
and introducing ablation is very small. The insensitivity of the static 
stability to gas composition is also demonstrated. The stability decreases 
with increasing angle of attack which indicates that the variation of pitching 
moment with angle of attack is not linear for this configuration. 

Estimates of C m<x at zero angle of attack calculated by means of conical- 
flow theory (ref. 11) and Newtonian theory (ref. 12) are also included in this 
figure. Comparison of these theoretical estimates (valid only at zero angle 
of attack) with the experimental data shows both theories underestimate the 
initial stability by about 20 percent. 
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Several methods are available (refs. 13 to 17 ) 'which permit analysis of 
the observed pitching and yawing motions of a symmetrical body having a non- 
linear pitching moment to obtain C m as a function of a. The method of 
reference 17; developed under the assumption that the nonlinear moment can be 
described by an arbitrary integer -power series of the resultant angle of 
attack, was used to analyze the data obtained in the present tests. Applica- 
tion of this method to the data of this report is discussed in appendix A. 

The local variations of pitching moment with angle of attack, extracted 
from the experimental data shown in figure 8 by the above procedure, are pre- 
sented in figure 9 * Examination of the experimental data at the three test 
Mach numbers shows, as expected, that the stability of the configuration is 
little affected by changes in Mach number from 4.5 to 13 .5 • The small 
increase in stability indicated at the higher angles of attack for an increase 
in Mach number from 4.5 to 9 .0 could very well be considered within the accu- 
racy of the experimental data at these angles of attack. The variation of 
pitching moment with angle of attack given by Newtonian theory and the initial 
slope calculated using conical -flow theory are also included in figure 9- As 
mentioned earlier, both theories underestimate the pitching moments at the low 
angles of attack by about 20 percent; however, Newtonian theory predicts an 
almost linear variation of C m with a; whereas the experimental variation is 
seen to be slightly nonlinear with angle of attack. The agreement between 
experiment and theory, therefore, improves with increasing angle of attack and 
is considered close throughout the entire angle -of -attack range of this 
inve st igat ion . 


Dynamic -Stability Characteristics 

The results of the dynamic -stability measurements are presented in 
figure 10, where values of the dynamic -stability parameter, |, are plotted as 
a function of the maximum resultant angle of attack, cr m . Negative values of 
5 represent a convergent model motion (dynamic stability) . These data show 
that the configuration is dynamically stable throughout the angle -of -attack 
range investigated for each Mach number and that it exhibits the same dynamic 
stability in air as it does in the carbon -dioxide --nitrogen mixture. These 
data also show that, within the scatter of the data, the values of | remain 
constant with increasing angle of attack and increasing Mach number . A value 
of | of -2.0 is a reasonable average for most of the data in figure 10, and 
is approximately equivalent to a convergence of about 7 percent per cycle . 
Examination of the data in figure 10 shows an appreciable amount of scatter in 
the values of | obtained at a Mach number of 13*5 for angles of attack below 
10°. At these conditions, three flights gave values of £ close to - 7 - 0 , 
which is approximately equivalent to a convergence of about 22 percent per 
cycle, and one flight gave a value of £ of about +2.0, which is equivalent 
to a divergence of about 7 percent per cycle . To determine the largest rea- 
sonable variation in the experimental values of £ for these flights, a 
probable error of 0.25° in angle -of -attack measurements, determined statisti- 
cally from many readings by several observers, was introduced to the input 
data of these flights. These errors in the worst possible arrangements 
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changed the values of £ by less than ±1.0. Therefore, these values of £ 
are considered representative of the actual damping characteristics experi- 
enced by the models in these flights. 

There are several factors which might have contributed to the scatter 
of these data. One such factor is of the possibility of model damage. Calcu- 
lations indicate that the stress at the bimetallic joint of the models could 
have, under an unusually high acceleration launch, equalled the compressive 
strength. The shadowgraphs from all the flights at Mach number 13-5 under 
high magnification, showed no evidence of deformation at the bimetallic joint 
but a very slight discontinuity in the bow -shock -wave angle of some of the 
models just ahead of the cylinder indicated there were surface def ormations . 
Although a compressive failure of the titanium should not have caused this 
discontinuity since the stress here should have been less than that at the 
bimetallic joint, the sabots which grip this section of the models (see 
fig. 3(c)) could have caused small irregularities in the plastic coating. It 
was not possible, however, to correlate this effect with the values of £ 
which produced the large scatter. Another possible explanation is that the 
scatter is due to effects of nonuniform ablation over the surface of some of 
the models and/or variations in the amount of ablation from flight to flight. 
Mass addition into the flow field about the body which resulted from ablation 
of the plastic is visible in the wakes in figures 5(c) and 6(c) . It has been 
shown In reference 18 that mass addition could cause large changes In the 
static stability, and consequently in the dynamic stability, of a body moving 
at hypersonic speeds and that these changes in stability Increase with 
increasing mass addition and depend on the angle of attack of the body. For 
the present tests, it was inferred from the experimental data presented in 
figures 7 and 8 that ablation effects were not sufficiently large to signifi- 
cantly affect the static stability and drag characteristics of the configura- 
tion; however, the dynamic stability might have been influenced. One other 
possibility is that the l-l/2 cycles of oscillation at Mach number 13*5 (see 
fig. ^4-) constitute a minimal amount of motion for defining the dynamic - 
stability parameter when o m is small, especially when "trim” is present. 

Some or all of these possibilities could have contributed to the scatter in 
the dynamic -stability results at this Mach number. 

Values of the damping -in -pitch derivative, C m ^ + C m( ^, were calculated by 
means of equation (5) using the values of £ (fig* 10) and Cp (fig. 7) meas- 
ured from the present tests, with values of CL a estimated by conical-flow 
theory, and are presented in figure 11 as a function of cr m . These data show 
the same trends as those observed for the values of £ in figure 10, namely, 
constant dynamic stability for the configuration throughout the Mach number 
and angle -of -attack range investigated. It is Interesting to note that the 
values of Cm. + C m ^ have the same sign as the corresponding values of £ 
and are about one -tenth as large. Thus, as seen from equation (5) the combi- 
nation of a stabilizing lift -curve slope (positive value of CL a ) and the 
stabilizing damping -in -pitch derivative overshadow the destabilizing contri- 
bution of drag to produce a convergent motion in constant speed flight at 
constant altitude . 
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Stability Characteristics for Large Amplitudes of Oscillation 


The present configuration is intended to be statically stable about only 
one trim attitude, nose forward about zero angle of attack. Although experi- 
mental measurement of the amount of static instability (positive value of 
Cm a ) provided for in the design of the configuration about the base -forward 
trim attitude of l80° was beyond the scope of the present investigation, two 
exploratory flights were made with the models launched backward to determine 
whether they would remain flying backward or, as expected from calculations, 
would begin righting themselves to a nose -forward attitude. The models were 
launched at an initial angle of attack near l80° (see fig. 3(d)) at the same 
test conditions mentioned previously for the other flights at Mach number 9 - 0 , 
namely, nominal model velocity of 3 .05 km/sec and test-section static pressure 
of 4 X 10 4 newton/m 2 ( 300 mm of mercury) . 

The measured time history of the angular orientation for one of these 
flights (flight 703), presented in figure 12, shows that although there was 
only about one -half cycle of oscillation in the available 12 .2 -meter -length 
test section, the model did not remain in the base -forward orientation and did 
not tumble but achieved a nose -forward orientation and appeared to be oscil- 
lating about zero angle of attack with amplitudes of oscillation as large as 
165°. This result indicates that the configuration is statically stable only 
about the nose -forward trim attitude of zero angle of attack. Another inter- 
esting and important possibility is suggested by this result and by the knowl- 
edge that for flight at constant or slightly decreasing dynamic pressure a 
model having negative or neutral dynamic stability would be expected to 
tumble . Since the model in this flight did not tumble , it is indicated that 
the dynamic -stability characteristics of the configuration, over this large - 
amplitude range, are such as to produce a convergent motion. 

To investigate the dynamic -stability characteristics of the configuration 
over this large -amplitude range and to substantiate the result of flight 703 
that the configuration has a unique nose -forward stable attitude, the second 
flight (flight 663) was conducted in the Ames Pressurized Ballistic Range 
(in air) which enabled the trajectory of the model to be recorded over a 
62-meter length in 24 shadowgraph stations, thus providing several cycles of 
oscillation. The measured time history of the angular orientation for this 
flight presented in figure 13* along with that of flight 703* shows, as did 
flight 703, that the model did not remain in the base -forward orientation and 
did not tumble but again achieved a nose -forward orientation and oscillated 
about zero angle of attack with large amplitudes of oscillation, (it should 
be noted that the model experienced about a 68-percent decrease in dynamic 
pressure during this flight.) Further examination of the experimental data in 
this figure shows that although the time history for flight 663 is not as well 
defined as that for flight 703, 1 the model in this flight achieved almost 
exactly the same peak amplitude at the first peak of oscillation as did the 


^■Although the wave length of oscillation was approximately the same for 
both flights the observation stations in the PHFF facility are spaced at 
1.22 -meter intervals while the shortest spacing between the observation 
stations in the PBR is 2.13 meters. 
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model in flight 703 • This almost exact duplication of the motions for both 
flights provides strong experimental evidence that the configuration is stati- 
cally unstable about the trim angle of attack of 180 and -will begin to over- 
turn, immediately, to a nose -forward orientation and oscillate about its 
unique stable trim attitude of zero angle of attack. Since evaluation of the 
stability derivatives from analysis of the high -amplitude experimental data of 
flight 663 could not be accomplished using the method presented earlier 
(ref. 8) to extract these derivatives from low -amplitude motions, an attempt 
was made to generate a theoretical motion which would provide a good fit to 
the experimental data and allow the desired stability derivatives to be 
inferred. The differential equation of motion of a missile oscillating in 
pitch at constant altitude (see ref. 10), written with respect to time as 

a - ^ ka.a - qMC m = 0 ' (7) 

where 

fa ■ + (s + vCC* (s> 

was programmed for numerical integration by a computer. Equation (7) "was 
iterated upon to obtain the best match to the experimental data. The time 
variation of pitching moment with angle of attack given by Newtonian theory, 
the time variations of dynamic pressure and velocity measured for this flight, 
along with various values of k x were used as inputs for these calculations. 
These calculations were initiated at the first peak value of a, l6^°, where 
the value of d is zero. The theoretical motion obtained in this manner 
which best fit the experimental data is presented in figure 13 and is seen to 
be an excellent representation of the actual motion experienced by the model 
during this flight. This excellent agreement between the theoretical motion 
and the experimental data indicates that the pitching moments predicted by 
Newtonian theory for this configuration are very close to the actual pitching 
moments experienced by the model throughout the angle -of -attack range from 0° 
to l80° . ( It should be noted that according to Newtonian theory the conf igu- 

ration is stable only in the nose -forward attitude.) Further examination of 
the data presented in figure 13 shows that the theoretical motion is conver- 
gent (the value of ki was - 1.15 and, as can be seen, is equivalent to a 
convergence of about b percent per cycle), which indicates that the configu- 
ration is dynamically stable, about zero angle of attack, throughout the 
entire angle -of -attack range from 0° to 180° . 


Application of the Present Results to a Full-Scale Vehicle 
Passively Entering the Atmosphere of Venus 

To determine the significance of the present results when applied to a 
full-scale vehicle passively entering a planetary atmosphere, calculations 
were made of the oscillatory behavior of the present configuration, initially 
oriented base forward, flying an example entry trajectory through a model of 
the Venusian atmosphere . Studies of the oscillatory behavior of missiles 
entering the atmosphere on ballistic trajectories (ref. 19 ) have shown that 
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the rapid increase in atmospheric density experienced by such vehicles 
strongly damps their oscillations. For the present example, two entry trajec- 
tories were considered: one, using an initial flight -path angle, of 89° 

which corresponds to a steep entry; and the other using a 7^ of 40° . These 
values for and the following assumed initial conditions were used to 

compute the entry trajectories for this example: 

Vj_ = 11. ^3 km/sec 
y^ = 160 km 
m/CpA = 156.8 kg/m 2 

The atmosphere used in these estimates was one selected on the basis of infor- 
mation presented in reference 20. As shown in reference 21, this atmosphere 
can be closely approximated by the well-known assumption that the atmospheric 
density varies exponentially with altitude 


for the altitudes between l6o km and 40 km (the altitudes of interest in the 
present example). This assumption and the following corresponding values of 
these parameters for Venus (see ref. 21) were applied to the present example: 

Py = I.58XIO 4 per m 

P Q = 531 kg/m 3 

These conditions define the entry trajectories considered. Some of the tra- 
jectory parameters, namely, the variations of velocity and dynamic pressure 
with altitude, were computed using the above conditions in equations derived 
in reference 22 and are presented in figure ±k . It can be seen from this 
figure that maximum dynamic pressure occurs at an altitude close to 65 km for 
both trajectories. 

Analyses which describe the oscillatory motions developed by vehicles as 
they descend through a planetary atmosphere are presented in references 23 
and 2 k . The analysis presented in reference 23 allows for an arbitrary ini- 
tial angle of attack (nose forward or base forward) but assumes the aerody- 
namic damping of the vehicle to be zero. The analysis presented in 
reference 2^ allows for aerodynamic damping of the vehicle but assumes con- 
stant aerodynamic coefficients which requires the vehicle to be initially 
oriented nose forward at a low angle of attack. Since the vehicle, for the 
present example, is to enter the Venusian atmosphere initially oriented base 
forward, the analysis of reference 23 was used to describe the oscillatory 
motions for the portion of the trajectory for which the aerodynamic damping of 
the vehicle could, in fact, be considered negligible compared to damping due 
to increasing dynamic pressure (see ref. 19) . This analysis showed the vehi- 
cle to be oriented nose forward and oscillating at low angles of attack at a 
point in the trajectory prior to reaching q max . ^ vas therefore possible to 
use reference 2 k to determine the effect of the measured damping on vehicle 
motions for the remainder of the trajectory. 
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Reference 23 shows that under the assumptions made in the analysis the 
envelope of these oscillations, expressed in terms of significant vehicle and 
planetary properties, can he written as 


(tan e)env _ qq 
tan 6p 

For the present examples, the vehicle was assumed to he initially oriented 
hase forward at an angle of attack of 179°.? with zero initial pitch rate, 
which gives values of the parameters C and G of 0.79^ and 0.0215, respec- 
tively. (The parameters C and G are dependent on the initial angle of 
attack and initial pitch rate assumed and cannot he evaluated explicitly; 
their values are determined from charts presented in reference 23.) The 
following physical properties of the vehicle were assumed 

A = 0.771 m2 

d - 1.0 m 

Iy = 6 -kk kg-m 2 

m = 6 .66 kg 

C D = 0.55 

Snax “ "°- 08 35 


-Pv- 2 sin 2 7-jIy Pyh (C D p 0 A/fVn sin 7i) e 


-!Vh 


1/4 


2p MC m 
*o %iax 


(9) 


where C m;max is the maximum value of pitching moment estimated using 
Newtonian theory. Values of the envelopes of oscillations, (a) e nvj computed 
are presented in figure 15 as a function of altitude . The results in this 
figure show that for that portion of the entry prior to q max the motions, 
for both entries, converge very rapidly so that at the altitudes where heating 
is most severe (between 80 km and 65 km, see fig. 1^-) , the amplitudes of 
oscillation of the vehicle are less than 35° • These results also show that 
for the altitudes subsequent to q max , an analysis which assumes the aerody- 
namic damping of the vehicle is zero predicts highly divergent motions. How- 
ever, as was shown in figure 10, the present configuration possesses positive 
aerodynamic damping and the effect of this positive damping on the oscillatory 
motions must he included in the analysis. 

Reference 2b shows that for constant aerodynamic derivatives the envelope 
of oscillations has the simple form 

n -K/2 

Cm z . 

°teax = , /4 

^00 
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It also has been shown in this reference that the dynamic -stability parameter 
K, in the form 




+ (C m(1 


Cm a,) 



( 11 ) 


is a convenient parameter for describing the dynamic stability of a vehicle 
descending through a planetary atmosphere. It can be seen from a comparison of 
equations (5) and (ll) that K can be expressed in terms of | and Cp as 
follows 


K = 


JL 

CD 


1 


( 12 ) 


Values of K computed from the measured values of £ (fig- 10) and Cp 
(fig. 7) are presented in figure l 6 as a function of a m . From these data it 
was considered that a value of K of -3-0 would be a conservative mean value 
to use in equation (l0). This value of K and the example entry trajectories 
presented earlier were introduced into equation ( 10 ), expressed in ratio form, 
and values of the amplitude ratio, (a/oci) max , were computed. These calcula- 
tions were initiated at an altitude of 80 km, since, as can be seen from fig- 
ure 15 , the oscillatory angles of attack of the vehicle below this altitude 
are less than 36 ° and, as demonstrated in figures 9 and l 6 , the aerodynamic 
derivatives of the configuration can be considered reasonably constant in this 
angle -of -attack range. Values of (a) env were computed from the ratio 
(a/ai) max using the angles of attack at 80 km given by the analysis of refer- 
ence 23 as initial angles of attack ( 36 ° for the 7± = 89° case and 30 ° for 
the 7 ^ = 40° case) and are presented in figure 15- These results show that at 
altitudes below that for q max ^ e effect of the measured aerodynamic damping 
of the configuration overshadows the effect of decreasing dynamic pressure to 
produce a strongly convergent motion. The amplitudes of oscillation at an 
altitude of 56 km for both examples are seen to be less than 5° • Therefore, 
from these calculations it is indicated that the present configuration, ini- 
tially oriented nose forward or base forward, should experience highly conver- 
gent pitching oscillations during entry into the postulated atmosphere of 
Venus . 


SUMMARY OF RESULTS 


The results of an experimental investigation to determine the stability 
and drag characteristics of a conical configuration at Mach numbers of 4.5, 
9-0, and 13.5 and Reynolds numbers from 0.4 to 1.2 million in a N 2 -C0 2 mix- 
ture and in air are as follows: 

1. The static stability, dynamic stability, and drag of the configura- 
tion are the same in air as in a gas mixture composed of 9 “Percent carbon 
dioxide and 91 - percent nitrogen for the entire range of Mach numbers and 
angles of attack of this investigation. 


17 



2. The configuration is statically stable about an angle of attack of 0° 
(nose-for-ward orientation) for the assumed location of center of gravity; 

0.72 d from the nose. Strong evidence is presented -which shows that the 
configuration is statically unstable about a trim angle of attack of l80° 
(base-forward orientation); so that; if initially oriented base forward; it 
will reorient itself to a nose-forward attitude. The static stability is 
essentially invariant with Mach number and slightly nonlinear with angle of 
attack in the angle range up to 40° . 

3 . The nonlinear variation of pitching moment with angle of attack was 
closely approximated by a two-term power series of the resultant angle of 
attack . 

4. The configuration is dynamically stable for steady flight at constant 
altitude. For this condition and for initial amplitudes of oscillation as 
high as 40°; the models experienced pitching motions which converged at the 
rate of about 7 percent per cycle. Also; there is strong indication that the 
configuration is dynamically stable at pitching amplitudes approaching 180° so 
that it will oscillate about the nose -forward stable attitude with diminishing 
amplitude of oscillation. Analysis of the large amplitude oscillations of a 
model concurrently experiencing a rapid decrease in dynamic pressure with time 
showed a convergence of about 4 percent per cycle. 

5. The drag coefficient of the configuration decreases approximately 8 
percent for an increase in Mach number from 4.5 to 9*0 but remains constant 
for a further increase in Mach number to 13*5* The drag increases approxi- 
mately 30 percent; for each Mach number; as the angle of attack is increased 
from 0° to about 28° . 

6. The variation of pitching moment with angle of attack given by 
Newtonian theory is in close agreement with the measured variation for the 
angle -of -attack range from 0° to about 36 °; although Newtonian theory and 
conical -flow theory underestimate the initial slope by about 20 percent . 
Furthermore; it is indicated that Newtonian theory gives close estimates of 
the pitching moments of this configuration for the entire angle -of -attack 
range from 0° to 180° . The drag of the configuration was accurately estimated 
by conical-flow theory at zero angle of attack and was closely approximated by 
Newtonian theory for the angle -of -attack range investigated; 0° to about 28°. 

7 . Calculations indicate that the measured stability characteristics 
of the present conf iguration should be adequate for orienting the vehicle 
properly during a Venus entry. For any initial orientation of the vehicle 
(nose forward or base forward) the pitching motions should converge to a very 
small fraction of the amplitude at entry. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field; Calif.; Dec. 18; 1964 
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APPENDIX A 


APPLICATION OF NONLINEAR PITCHING -MOMENT ANALYSIS 


A method which permits analysis of the pitching and yawing motions of a 
spinning symmetrical body having a nonlinear pitching moment to obtain the 
pitching moment as a function of angle of attack is described in reference 17* 
The method allows for a more general reduction of free -flight data by assuming 
that the restoring moment can be described by an arbitrary power series of the 
resultant angle of attack, of the form 

(Al) 

n=o 


M( a) 


CO 

- - £ 


where 


M(a) = ^ Cjn 


(A 2 ) 


A simple expression is then derived which relates the frequency of oscillation, 
w, obtained from data analysis with linear equations, to the maximum and mini- 


mum amplitudes for an infinite number of possible moment combinations, 
expression, which includes the nonlinear inertial terms due to spin, is 
ref . .17) 
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By use of this expression, values of the re storing -moment coefficients, M n , 
can be obtained from the observed frequencies of oscillation and amplitudes of 
independent flights of the models at similar flight conditions but with differ- 
ent amplitudes of oscillation. Equation (A3), as discussed in reference 17 , 
is an approximate solution which becomes exact in the limit of circular motion, 
that is, equal pitch and yaw, 90° out of phase and is least accurate for the 
case of planar motion. As can be seen from examination of figure 4 , the 
motions to be analyzed for the present investigation are fairly planar. How- 
ever, the applicability of the method to the case of planar motion was 
examined in reference 17 for certain assumed moment representations for which 
exact solutions are possible . Comparison showed that the approximate solution 
(eq. (A3)) gives results extremely close to the results obtained from exact 
solutions. Further comparison in reference 17 shows, again for the case of 
planar motion, that the approximate solution gives results that are, for the 
most part, more accurate than results obtained by the technique of refer- 
ence 13. Therefore, a high degree of confidence in the applicability of this 
method to the data of this report is justified. 
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Application of the method to the present experimental data *was 
accomplished in the following manner. Equation (A3), which allows for an 
infinite number of possible moment combinations, was limited, from practical 
considerations, to the first seven terms written as 


5 5 

°m + °o 


^ = Mo + | M x 


°m + °( 
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(°m 2 - °o s ) ‘ 
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32 


(A 4 ) 


(For convenience, since values of are available from the linear analysis 

of the experimental data, the substitution w 2 = (qAd/ly)C ma was made in 
equation (A 4 ).) An examination of the experimental data in figure 8 indicated 
that the nonlinearities could not be considered large; therefore, it was 
decided that a two- or three -term moment in the resultant angle of attack 
would adequately represent the moment curve. Equation (Ah) was programmed for 
machine computation so that all possible moment combinations, each containing 
the linear term, M Q , plus one and two members of the set (Mi, M 2 , M3, M4, M 5 , 
Me) were fitted to the data by the method of least squares, The results are 
presented in figure 17, which shows the pitching -moment curves given by the 
various assumed moment representations that best fit the experimental data for 
each Mach number. The linear plus cubic pitching -moment curves are also 
included in this figure. (These pitching -moment curves are the same as would 
be obtained by the method of reference 15, which develops this case exclu- 
sively.) It is evident from these results that for the angle -of -attack range 
investigated the moment curves given by the best two- and three -term moment 
approximations for a particular Mach number are close enough together that it 
is unnecessary to choose between them. Furthermore, since the results show 
that the nonlinear pitching moment can be closely approximated by the simplest 
of the moment representations tried, namely , a two-term moment representation, 
the applicability of the method to the data of this report rests on the 
approximate linear relationship of C m<x with the maximum and minimum ampli- 
tudes of the motion. Therefore, in order to determine the validity of this 
assumption, values of C ma were plotted as a function of the appropriate 
amplitude term given in equation (A4-) for a one-three moment representation at 
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Mach member 4.5 and for a one-four moment representation at Mach numbers 9*0 
and I 3.5 for each flight and are presented in figure 18. It is evident from 
this figure that the data fall very closely along straight lines for each 
Mach number and, therefore, the particular two-term moment representations of 
the method are excellent approximations. 
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(c) Nominal Mach number = 13*5 
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(a) Model. 



(t>) Plastic-coated model. 

Figure 1.- Sketches of models showing nominal dimensions. 


First station. 


(a) N 2 -C0 2 


Last station. 



First station. Last station. 

(b) Air 


Figure 2.- Shadowgraphs showing negligible tip blunting on ablating plastic-coated models; M 
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(c) Flight 8o6; M = 9.0 (d) Flight 809; M = 9-0 

Figure k.- Continued. 
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(e) Flight 987 ; M = 13-5 

Figure 4.- 



(f) Flight 985 ; M = 13.5 
Concluded . 
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(c) M = 13.5 

Figure 5*- Shadowgraphs showing flow-field phenomena of models flying in a 
9-percent carbon dioxide, 91 ~P e ^cent nitrogen mixture. 
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Figure 6.- Shadowgraphs showing flow-field phenomena of models flying in air. 
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Conical-flow theory, M=4.5 
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Newtonian theory 


Conical -flow theory, M = 9.0 and 13.5 
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Solid symbols indicate air data 


Root -mean -square resultant angle of attack, cr 

Figure 7*- Drag results. 
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Figure 8.- Static stability results (linear analysis of data) , x C g/d = 0.72. 
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Figure 9.- Comparison of experimental pitching-moment coefficient with theory; x C g/d = O.72 
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Figure 10.- Effects of Mach, number and test gas on the dynamic stability parameter, %. 
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Figure 11.- Effects of Mach number and test gas on the damping -in-pitch derivative, C 
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Resultant angle of attack, 
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Figure 13 -- Fit of theoretical motion equation (eq. (7)) to experimental data; x c g/^ - 0 .72 
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Figure 15*- Os dilation -amplitude histories of the configuration entering a 

model Venusian atmosphere. 
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Figure 1 6 .- Effects of Mach number and test gas on the dynamic stability parameter, K. 












Best 2-term moment approximation (1-3) 


Best 3-term moment approximation (1-5-6) 


Angle of attack, a, deg 


(a) M = 4.5 

Figure 17 .-* Pitching-moment coefficient corresponding to the "best 2- and 3-term assumed moment 

representations j x C g/d = O .72 d. 




(1-3) moment approximation 
Best 2-term moment approximation (1-4) 


Best 3-term moment approximation ( 1-2-3) 
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Figure 17*- Continued. 
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Best 2-term moment approximation (I -4)- 
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"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 
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